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DEVELOPMENT OF A SIMULATOR FOR STUDYING 
SIMPLIFIED LUNAR ESCAPE SYSTEMS 

By George J. Hurt, Jr., David B. Middleton, 
and Marion A. Wise 
Langley Research Center 

SUMMARY 

A lunar-escape -system simulator was designed and fabricated to be used as a 
research tool to study problems related to lunar escape and simplified manual guidance 
and control. The primary requirements for the simulator were based on items that 
would be desirable for an actual lunar escape system: first, simplified propulsion; 
second, a simplified stabilization and control method; and third, a simplified guidance 
and navigation technique. To ensure exact trajectory calculations the real-time digital 
computer program used in the simulation was compared with an independent computer 
program, and good agreement was obtained. 

The facility design included load cells to transmit pilot kinesthetic control motions 
realistically. Most of the vehicle characteristics and control logic were contained in the 
computer program; thus investigation of a wide range of vehicles (including the lunar 
flying vehicles) and guidance and control systems was possible. During checkout opera- 
tions test subjects were able to adapt to the equipment with a minimum of training and 
instruction in the operation of the simulator. 

INTRODUCTION 

The ascent stage of the lunar module (LM) includes many redundancies and other 
features to maximize its chance of successful operation. It must, however, rely on a 
single rocket engine for lift-off; an alternate means of escaping to lunar orbit is thus 
highly desirable. A number of approaches have been suggested including multipurpose 
surface -to-orbit shuttlecraft and lunar flying vehicles (LFV) with the capability of also 
achieving lunar orbit. Both of these approaches involve rather extensive development 
time and probably a separate earth launch to place the vehicles on the moon. 

A third approach which has for several years been considered at NASA Langley 
Research Center (LRC) is that of packaging a very simplified lunar escape system (LES) 



on the LM (see refs. 1 to 4). Primary emphasis has been on the development of appro- 
priate guidance and control techniques and determination of stability requirements for a 
simplified LES. Some effort at LRC has also been directed toward development of a 
dual-purpose LFV-LES. The problem of unloading, setting up, and servicing an LES 
(or LFV) on the moon is also discussed in these references. 

Flight-simulation studies are now essential in determining the feasibility of using 
particular minimum -complexity LES configurations, control modes, and guidance tech- 
niques. A lunar-escape-system simulator (LESS) has been developed at LRC to study 
such problems. Throughout this report the lunar escape system will be referred to as 
the LES; the simulator of that system will be LESS. The primary requirements for the 
LESS were based on items that would be desirable for an actual lunar escape system: 
first, simplified propulsion; second, a simplified stabilization and control method; and 
third, a simplified guidance and navigation technique. 

The LESS is the subject of the present report, which covers the development of 
simulation hardware and an associated real-time digital computer program. Most of the 
vehicle characteristics and control logic are contained in the computer program; thus 
investigation of a wide range of vehicles (including the LFV) and guidance and control sys- 
tems is possible. The real-time computer program has been compared with an inde- 
pendent digital computer program with good agreement. The LESS is operational at 
present and has proved to be a useful representation of the pilot problems involved in 
using either jets or kinesthetic inputs for attitude control. 

SYMBOLS 

Values are given in both SI and U.S. Customary Units. The measurements and 
calculations were made in U.S. Customary Units. 

A e cross-sectional area of exit nozzle of rocket main thruster 

a semimajor axis of vehicle orbit 

a^j direction-cosine elements of matrix used for transformation from body 

to inertial system (i=l ,2 ,3; j=l ,2 ,3) 

b^j direction-cosine elements of matrix used for transformation from body to 

local-vertical system (i=l ,2 ,3; j=l ,2,3) 

Cjj direction-cosine elements of matrix used for transformation from inertial 

to local-vertical system (i=l,2,3; j=l,2,3) 
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Dj,D 2 ,D 3 auxiliary variables used to simplify moment equations in appendix B 
F magnitude of F (or gage thrust) 

Fj,F 2 maximum and intermediate constant magnitudes of F 
Fft magnitude of F^ 

F XB’ F y B ’ F z B body-axis components of F 

Fx B ’Fy B >Fz B body-axis components of F 

F force due to rocket engine (or to main thruster of LES, in particular) 

F gxt external force acting on rocket mass 

F^ horizontal restraining force of rocket test stand 

F b total restraining force of rocket test stand 

F force due to combination of thrust and gravity 

G vehicle reference point (viz., initial location of vehicle center of gravity (c.g.)) 

g e acceleration due to earth gravity, 9.81 m/sec^ (32.2 ft/sec^) 

g m acceleration due to lunar gravity, 1.62 m/sec^ (5.32 ft/sec^) 


h a ,hp altitude of apocynthion and pericynthion, respectively 

I s p specific impulse of LES propellant mixture 

IxxJyyJzz moments of inertia of vehicle with respect to body axes 


I X y,I X z>Iyz products of inertia of vehicle with respect to body axes 



moments of inertia of mj with respect to temporary axes X, Y, 
and Z through the c.g.' s of mj (j=l,2,...,5) 
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products of inertia of mj with respect to temporary axes X and Z 
through the c.g.’s of mj (j=l,2,...,5) 


( dIxx )disk , ( dIzz )disk 


differential moment of inertia of elemental propellant disk with 
respect to temporary axes X and Z through center of the 
disk 


m instantaneous mass of vehicle 

m i elemental mass particles which comprise m (i=l, 2,. ..,«>) 

mj mass modules which comprise m (j=l,2,...,5) 

mp- initial propellant mass 

0 jO 

Nq angular momentum of vehicle about its instantaneous c.g. (point O) 

N x B ’ N y B ’ N z B body-axis components of Nq 

O location of instantaneous c.g. of vehicle (also origin of body axes X B , 

Y B , and Z B ) 


P,q,r 


magnitude of p g 

directional rocket-exhaust pressure at nozzle exit 

components of angular velocity measured about body axes X B , Y B , and Z B , 
respectively 


Qt ^ »Qt tr .Q t rr body-axis torques due to attitude jets 

j j x B J ’ Z B 

Qp ,Qp miscellaneous body-axis torques 

-^> X B r,y B 

Q x ,Qyg,Q z body-axis components of total body torque 

R g ,Rq position vectors of point G and point O with respect to center of moon 
position vectors of m^ with respect to center of moon (i=l,2,...,°°) 
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Ro 


magnitude of Rq 


^B ,R yB ,Rz B 


body-axis components of Rq 
inertial components of Rq 


GO 


•m 


t l» t 2’ t 3 


l BO 


dV 


V z 


B 


0,h 


position vector of point O with respect to point G 

position vectors of m^ with respect to point G (i=l,2,...,°°) 

radius of moon 

radius of pericynthion 

radius of spherical propellant tank 

time 

time of beginning, midpoint, and end, respectively, of the 90° LES 
pitch-transition maneuver (pitchover) 

time of burnout 

body-axis components of Vq 

body-axis components of Vq 

elemental volume of disk of propellant 

simulated velocity-along -the -thrust-axis signal 

magnitude of v e 

horizontal component of Vq 

velocity of general rocket mass used in derivation of force equations 
effective exhaust velocity of rocket-engine gas particles with respect to m 
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v G total linear velocity of point G with respect to center of moon 

Vq total linear velocity of point O with respect to center of moon 

Vqo velocity of point O with respect to point G 

x,y,z distances (or directions) in X,Y,Z system 


x B’ y B’ z B distances (or directions) in Xg,Yg,Zg system 

x r .,y r -,z r . body-axis components of initial positions of p m . (j=l,2,...,5) 
^>3 j 

xpyp z I distances in XpYpZj system 

x i’ y i’ z i body-axis components of r. (i=l,2,...,°°) 

x LV’ y LV’ z LV distances in X LV ,Y LV ,Z LV system 

Z B t displacement distance of thruster gimbal below point G (measured along 

Zg-axis) 

z 1 distance from center of propellant tank to propellant surface 


out-of-plane angle 


Ax,Ay,Az body-axis components of r^ 

<5xg 3>®yg 3 body-axis components of control pilot’s c.g. shift 

6 Z 5 c.g. shift of propellant mass with respect to center of spherical tank 
B’ 

T) down- range angle 

(p,ip,d Euler angles used in transformation from body to inertial system 

cp T ,i// ,6' Euler angles used in transformation from body to local-vertical system 

6^,62 constant pitch angles used before and after pitchover, respectively 
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M 


lunar gravitational constant, 4.9028 X 10*^ m 3 /sec 2 
(1.7314 x 10 14 ft 3 /sec 2 ) 


Pm j 


Pjr. »Py. >Pz. 
X ] y ] Z ] 


body-axis components of the thrust-misalinement angle 

position vectors from the vehicle instantaneous c.g. to the component 
masses mj (j=l,2,...,5) 

body-axis components of p m ^ (j=l,2,...,5) 

propellant mass density 


w B > x I » w B,y I > a, B > z I 


inertial-axis components of body-axis rotational rate 


w T xr „ ,co T ,co inertial-axis components of rotational rate of the local 

J-iV,xj LVyj’ LV,zj 

vertical 


components of o) measured in local-vertical axis system 

X LV VLV z LV 

co total rotational rate of vehicle about point O 

Axis systems (all systems initially parallel): 


X B ,Y B ,Z B body axes with origin at instantaneous c.g. of vehicle (point O) 

X G ,Y G ,Z G initial location of X B ,Y B ,Z B (point O s point G) 

Xj,Yj,Zj inertial axes (Zj lies along local vertical at lift-off site, and the XjZj-plane 

coincides with plane of command-service-module orbit) 

X B y,Y B y,Z B y local - ve rtical axes with origin at point O (Z B y is parallel to 

local vertical, and Y G v is parallel to Yj) 

X,Y,Z temporary axes with origin placed in convenient locations to aid in inter- 

mediate calculations (such as moments of inertia of mj (j=l,2,...,5) 
with respect to c.g. of each mj) 

A dot over a variable denotes differentiation with respect to time. A matrix 
inverse is denoted by £ ^J - ^, body-axis components of a vector are denoted by . 
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DESIGN CONSIDERATIONS 


In a real emergency lunar escape situation, it seems reasonable (see ref. 5) to 
require the astronauts to perform a number of guidance and control functions which 
would normally be handled automatically by LM onboard systems. This approach would 
lead to minimum-complexity hardware, which allows (1) possible stowage of the LES on 
the LM, (2) a relatively shorter LES development time, and (3) the high reliability asso- 
ciated with simplified systems. In particular, it seems desirable to limit the dry mass 
of the LES to approximately 135 kg (300 lbm); this limitation suggests that use must be 
made of the LM propellants (which become available in an emergency escape situation). 
The minimum weight of a suitable hypergolic propulsion system (engine, tanks, pres- 
surization system, etc.) could exceed one -half of this dry-mass allotment; this require- 
ment dictates a very simple and lightweight LES vehicle structure with very little 
instrumentation. 


Mission Requirements 

The basic mission requirements postulated for an emergency escape situation are 
the following: 

(1) The LES must lift two astronauts from the lunar surface to a safe lunar orbit. 
(By definition a safe lunar orbit must have a pericynthion altitude greater than 

15 240 meters (50 000 feet).) 

(2) The orbiting command service module (CSM) must rendezvous with the LES 
and pick up the two LES astronauts within the 4 -hour operating limit of their portable life 
support systems (PLSS). Another constraint is that the CSM must make the rendezvous 
within its energy allocation for LM rescue (at present, this allocation corresponds to a 
CSM characteristic velocity of approximately 240 m/sec (787 ft/sec). 

In addition to the basic mission requirements, the following assumptions are made 
for the LESS study: 

(1) The LES must be capable of being controlled by one astronaut. The second 
astronaut will be strictly control-inactive, although he may have some monitoring or 
communications duties. 

(2) On the lunar surface the LES will have c.g. sensing capability (e.g., load cells), 
which will allow the LES pilot to adjust the c.g, onto the expected thrust line for lift-off. 

(3) The required time of lift-off is known, either from communication with or visual 
sighting of the CSM. 

(4) Four propellant tanks are symmetrically arrayed with no propellant slosh 
(i.e., by use of positive expulsion bladders). 
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(5) The astronauts will use the LM life support provisions until just prior to LES 
lift-off. After lift-off the astronaut’s PLSS will be necessary in the open cockpit of 
the LES. 

An artist's concept of a LES is shown in figure 1. 



Figure 1.- An artist's concept of a LES. 

Guidance Considerations 

Guidance requirements of the LES are relaxed somewhat by the assumption that the 
CSM will perform the rendezvous. For this reason and also because of LES weight limi- 
tations, only unsophisticated guidance techniques are considered in the LESS studies. 

For example, holding the LES vehicle attitude to a series of constant pitch angles for 
specified time intervals is an open-loop method of launching into orbit, which requires 
only a clock and simple attitude indicators for instrumentation. An improved variation 
may be a three-axis attitude indicator (8 -ball) and the digital readout of an integrating 
accelerometer to provide velocity information. 

Control Techniques 

Several attitude control techniques requiring little or no control mechanization are 
being considered in the LESS studies. One is kinesthetic (balance -reflex) control of pitch 
and roll, which reduces attitude controller requirements from a three-axis hand controller 


L 



to a three -position on-off yaw switch (or stick). The type of kinesthetic control con- 
sidered for the LESS studies involves the ability of the control pilot to change and/or 
correct his attitudes by shifting his c.g. with respect to the vehicle line of thrust. From 
a standing position he can lean (or step) forward to shift his body c.g. (and consequently 
the total-system c.g.) forward of the thrust line and thereby induce a pitch-down vehicle 
torque. Similar torques are induced about the roll axis by leaning sideways (see fig. 2). 



Figure 2.- The kinesthetic motion used by a pilot 
to initiate pitch and roll moments with respect 
to the vehicle line of thrust . 

Kinesthetic control in various forms has been studied for a number of years (refs. 6 
to 8) and found to be possible but not desirable for most flight applications requiring 
precise control, particularly during landing. Lunar escape requires no landings; thus 
kinesthetic control may be feasible. In addition, it appears that kinesthetic control should 
require very little training because it makes use of the pilot's inherent ability to keep 
his balance. 

A second attitude control technique, thrust-vector control, requires very simple 
mechanization. A single main thruster can be double -gimbaled, and direct mechanical 
linkages between this thruster and a hand controller allow the pilot to tilt the thruster to 
generate pitch and roll moments. Jetavators or appropriate flaps on the thruster nozzle 
can be similarly linked to the controller for yaw control. With thrust-vector control, 
the pilot can either stand or be seated; if seated, he is less likely to introduce disturbing 
torques through inadvertent body shifts while operating the controllers. In some instances 
the pilot may wish to shift the total system c.g. intentionally in order to compensate for a 
bias torque (e.g., a torque due to thrust misalinement) . If standing, he can step in the 
appropriate direction; if seated, he can reposition his legs. 

The LESS is also appropriate for the study of other control techniques, such as 
using attitude jets only, or any combination of the techniques just mentioned. 
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LUNAR-ESCAPE -SYSTEM SIMULATOR 


General Description of LESS Hardware and Computer Program 

A lunar-escape-system simulator (LESS) was designed and fabricated to be used 
as a research tool to study problems related to lunar escape. The LESS was designed 
to accommodate a broad spectrum of studies related to lunar escape and simplified 
manual control. The primary requirements for the LESS were based on items that would 
be desirable for an actual lunar escape system: first, simplified propulsion; second, a 
simplified stabilization and control method; and third, a simplified guidance and naviga- 
tion technique. The basic simulator hardware consists of a fixed-base pilot control sta- 
tion appropriately tied into a digital computer which operates in real time. The fixed- 
base pilot control station includes provisions for (1) simplified thrust and attitude control, 
(2) instrument displays, and (3) equipment to sense the motions of a pilot's center of 
gravity (viz., strain-gage load cells). Features which can be added for particular studies 
are visual displays to the control pilot and free but limited pitch and roll motions of the 
control station itself. 

Figure 3 is a block diagram showing the man-hardware -software relationships of 
the LESS. This diagram illustrates that the LESS pilot control station can accommodate 
two men (standing or seated) but that it is outfitted primarily for one-man control. Con- 
trol inputs from the load cells and from the pilot's hand controllers are sent over tele- 
phone lines to the central computing complex where they enter the real-time computer 
through analog -to -digital signal converters. In turn, the computer updates the escape 
trajectory, calculates the instantaneous attitude of the vehicle with respect to a selected 
reference (e.g., the local vertical), and by means of digital -to -analog converters, returns 
selected analog signals to the simulator hardware. Several papers describing LRC 
simulation facilities are summarized in reference 9. One of these papers with direct 



Figure 3.- Block diagram of lunar escape system simulator (BESS). 
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application to the LESS is entitled "A Description of LRC Digital Computer Complex and 
Special Features for Real-Time Simulation Applications." The LESS was set up to oper- 
ate in the "multiprograming mode" described in this paper, and it "time-shared" with 
two other real-time simulations. A unique feature of the LESS system is some special 
circuitry and program logic which allows the pilot (at the remotely located pilot control 
station) to synchronize system operations by means of his thruster switch. For example, 
the computer can be used to calculate vehicle attitudes prior to lunar lift-off while the 
pilot is kinesthetically balancing the vehicle c.g. onto the expected line of thrust. Then, 
when thrust is turned on, the trajectory calculations begin and a discrete output signal 
from the computer starts the pilot's digital clock. 

More detailed discussions of the simulator hardware and computer program are 
given in the sections which follow. 


Pilot Control Station 

A simulated pilot control platform was clamped to the top of a sandwich-type load- 
cell structure. The relative size of the simulator assembly is illustrated in figure 4, 
which shows the pilot and passenger in position for a simulated lift-off. The control- 
platform equipment included a handrail, a yaw attitude controller, a main-engine thrust 
controller, and a situation display panel. 



L -69-4990 

Figure i+.- The simulator pilot-control station. 
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Figure 5 depicts the seated situation as provided for in the LESS design. A four- 
legged stool with a large footboard was placed on the LESS platform approximately over 
its c.g. For unbalanced pitching torques, the pilot can slide both feet forward or aft; for 
bias roll torques, he can shift one or both legs to one side. 



Figure 5 <- Sketch of pilot in seated position. 


Moment sensing .- A system of strain-gage load cells was used to sense the kines- 
thetic moments created by the pilot. The assembled load-cell sandwich and the amplifier 
rack are shown in figure 6. The orientation of the load cells is shown in figure 7. The 
main support post in the center of the platform acted as an XgYg-logic pivot and was 
designed to carry the static weight of the simulator assembly. The main support post 



L-69-1570.I 


Figure 6.- The assembled load-cell sandwich and the amplifier rack. 
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Figure 7.- A sketch of the load-cell orientation and associated hardware. 

incorporated a flexure section designed to permit sufficient hinging of the upper plate to 
allow maximum rated deflection of the load cells. Each of the load-cell connector links 
to the upper plate also included a flexure post to isolate the load cells from shear loads 
that might be generated by pitch and roll moment combinations imposed on the upper 
plate. A mechanical deflection stop on each load cell and a caging mechanism consisting 
of four corner posts and wedges were used to avoid inadvertent overloads during simula- 
tor alterations and shutdown periods. Large offcenter loads, such as a passenger or 
display packages, were counterbalanced with lead weights. Minor day-to-day out-of- 
balance moments were trimmed to zero via a trim adjustment in the load-cell amplifier 
system. 

Extraneous pickup from local 400-Hz and 60-Hz power lines was attenuated by 
filtering before transmission of the load-cell signals to the digital computer. This fil- 
tering also attenuated high-frequency components due to mechanical vibrations (the simu- 
lator was not mechanically isolated from the building). 

Thrust and attitude controllers .- The three-position LESS thrust control switch is 
shown in figure 8. The full-down switch position is thrust off; the center position and 
full-up position are wired to command various levels of intermediate and maximum 
thrust, respectively. 
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L- TO- 1344.1 

Figure 8 .- The simulator control platform. 


Two attitude controllers are provided at the pilot control station - a three-axis 
hand controller and a single-axis yaw stick controller. These controllers may be seen 
in figure 9. The three-axis controller shown is a CSM type (ref. 10). It is appropriate 
for most escape profiles in which attitude control is by thrust vectoring (gimbaled main 
engine) or in which there is all-jet attitude control. The yaw-output-signal conditioning 
was managed in the computer program. For the three -axis controller the attitude sig- 
nals could be simple on-off pulses or proportional to controller deflections. An electri- 
cal dead band was used in each of the three controller axes to ensure a zero electrical 
signal output for each mechanical control center. When vehicle pitch and roll moments 
were to be controlled kinesthetically, the pitch and roll potentiometer outputs of the three- 
axis controller were grounded to prevent spurious control signals from being introduced 
into the computer program. 

The single-axis yaw controller is a simple left-right spring -centered lever with 
microswitch contacts for signal initiation. The lever displacement required to generate 
a yaw input was fixed at a minimum. 
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It-' TO- 13^5.1 


Figure 9*- The three-axis hand controller and the single-axis yaw controller. 


Information displays .- In line with the simplified guidance and navigation concept, 
it was believed desirable to make some determination of the minimum attitude informa- 
tion necessary to enable the pilot to make a successful flight. The first information 
display group consisted of a digital clock, a fuel indicator, and a three-axis 8-ball. The 
instrument panel was hinged in order to accommodate the various eye -level heights 
which would be encountered (for the shortest subject seated to the tallest subject standing). 

Vehicle attitude information was presented to the pilot by a three-axis 8-ball. The 
instrument used was a prototype instrument from a lunar -module attitude information 
group (see ref. 10). A full-face presentation of the 8-ball is shown in figure 10. Various 
functions of the 8-ball assembly were either masked off or altered to suit the LESS study 
requirement. The roll steering needle was programed to present out-of-plane informa- 
tion to the pilot during the trajectory. A displacement of the needle indicated the direc- 
tion and relative magnitude of the out-of-plane velocity along the YB-axis. The pitch 
and yaw steering needles were not used for their normal function. They were coupled to 
the pitch and roll load-cell signals and used as vehicle leveling indicators for pitch and 
roll attitude during the prelaunch phase. The pitch and yaw needles were not used after 
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lift-off. Each of the three needles is shown in a displaced position in figure 10. A view 
of the instruments with the attitude rate information masked off is shown in figure 11. 



L-70-13V7.I 


Figure 10.- A pilot's view of the three-axis 8-hall in unmasked, configuration. 



L-70-1350.1 

Figure 11.- A pilot's view of the three-axis 8-hall with the attitude rate indicators masked. 
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LESS Computer Program 


The analog-digital real-time computing system was selected for the LESS studies 
because a high degree of resolution was necessary in computing the LES ascent trajec- 
tories and resulting orbits. The LESS computer program was written for six-degree- 
of -freedom motion referenced to the center of the moon. The control-input loop of the 
program had provisions for selecting several modes of attitude control individually or in 
combination. In particular, the kinesthetic control input signals are active for all 
research runs, so the pilot and passenger must stand (or sit) relatively still when thrust- 
vector control or attitude jets are being used; this same feature, however, allows the 
pilot to use kinesthetic control inputs to augment either of the other two control modes. 

The program was written in the FORTRAN IV computer language and had the option 
of selecting from four different numerical integration techniques. A sampling rate (and 
integration interval) of 1/32 of a second was used. This rate allowed sufficient accuracy 
for the trajectory calculations. Also, the digital outputs converted to smooth analog sig- 
nals for driving the pilot's displays (i.e . , the pilot saw 8-ball changes as continuous rather 
than stepped motions). 

Other elements of the computer program such as equations, reference axes, and 
axis transformations are discussed in the appendixes. Appendix A contains a description 
of the reference axes and equations for generation of the various axis transformations in 
terms of direction cosines. Also, Euler angle expressions were extracted from appro- 
priate elements of the transformation matrices; these expressions were then used to 
generate drive signals for the three-axis 8-ball. Appendix B contains derivations of the 
equations of motion in vector form and also in terms of body-axis components. The 
LESS inertia equations are constructed in appendix C, and finally equations for calcu- 
lating pertinent LES orbital parameters are developed in appendix D. 

Checkout of LESS Computer Program and Hardware 

The LESS received two types of verification tests - the computer program only 
and then a full system test with a pilot and all the hardware tied in. By using the equa- 
tions developed in the appendixes, the LESS computer program was test run in an ana- 
lytical mode by means of a special automatic control loop. 

To obtain initial conditions and other thrust and guidance parameter values for the 
test run, an independent digital computer program was written on the basis of the fol- 
lowing planar-motion equations: 
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Ro = Ro^ + |jCos0-^ 2 


"A 


RqV = ^ sin 6 - 2RqV 


m = - 


^spSe 


J 


( 1 ) 


where I S p is the specific impulse of the propellant and where the pitch angle 6 is 
referenced to the local vertical. The pitch angle was constrained as follows: 


0 = = 0 
_lo 

sec- 1 



(O^t^tj) 

(t x < t < t 2 ) 

(t 2 stst s)r 

( t 3 <t = * 30 ) 


( 2 ) 


The 90° pitch maneuver occurring between t^ and t 2 takes approximately 19 seconds. 
Also, the thrust level F is changed from to F 2 at t^ and remains at this 
second constant level until the time of burnout tgQ. 

To obtain a set of trajectory burnout conditions which would place the LES in a 
60-nautical-mile circular orbit, a series of iteration runs using equation sets (1) and (2) 
were made to determine values of Fj, F 2 , tj, and tgQ. For these runs an integra- 
tion subroutine involving variable integration steps based on specified limits for com- 
puter round-off and truncation errors was used. Parameter values obtained from this 
iteration process (planar program values) are shown in table I. These planar program 


TABLE I.- TRAJECTORY CHECKOUT PARAMETERS 


Parameter 

Values 

Automatic 
LESS program 

Planar program 

(F^ym) at t = 0, m/sec 2 (ft/sec2) . . . 

2.8032 (9.1968) 

2.8103 (9.22011 

(F 2 / 111 ) at t = tj, m/sec 2 (ft/sec 2 ). . . 

3.1504 (10.3360) 

3.1503 (10.3356) 

^/m) at t = tgQ, m/sec (ft/sec 2 ) . . 

5.4867 (18.0010) 

5.4864 (18.0000) 

tj, sec 

247.5 

246.6 

l BO> sec 

653.39 

652.54 
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values were inserted in the LESS program (with the automatic-control loop), which was 
run several times with different numerical integration techniques and fixed-step integra- 
tion intervals used. From these runs it was determined (1) that the automatic LESS 
program produced trajectory results in close agreement with those of the planar program, 
(2) that a fixed-step integration interval of 1/32 of a second was satisfactory, and (3) that 
the Runge-Kutta second-order (RKII) integration technique produced results almost identi- 
cal to those obtained by using the more sophisticated fourth-order Runge-Kutta method. 
Thus, to reduce computation time, the RKII method with an integration interval of 1/32 of 
a second was selected for use in the piloted LESS program. 

Before beginning the piloted LESS runs, the values in table I for the planar pro- 
gram were adjusted slightly to those shown for the automatic LESS program in order to 
obtain a more nearly circular reference orbit. The results of using both sets of values 
in their respective programs are shown in the tabular insert of figure 12. As indi- 
cated, neither solution produced an exactly circular orbit at a 60-nautical-mile altitude 
(60 n. mi. = 111 120 meters = 364 567 feet), but further iteration did not seem to be 
justified. 

The results of a typical piloted LESS trajectory are shown in figure 12. Time 
histories of the piloted LESS trajectory state variables are shown on the left side of this 
figure, and the corresponding burnout conditions and orbital altitudes are given in the 
last column of the tabular insert. The top four traces on the left side give attitude 
histories; the fourth trace down gives the time history of the pitch angle before the 90° 
pitch maneuver (pitchover), and the trace directly above it gives the time history of the 
pitch angle after the 90° pitch maneuver. Both pitch and roll were controlled kines- 
thetically, and yaw was controlled by means of small on-off jets. It can be noted that both 
pitch-angle and roll-angle excursions increased with time into the flight; these increases 
occurred because the vehicle became more responsive and harder to control as it lost 
the majority of its mass (the propellant mass constituted over 60 percent of the initial 
mass). 

Similar traces for the LESS operating with the automatic control loop are shown 
on the right side of figure 12. The roll and yaw angles were constrained to zero (not 
shown), and the pitch-angle history is shown on a single trace. Burnout conditions and 
orbital altitudes are given in the column of the tabular insert labeled "Automatic LESS;” 
these values were used as reference values for all the piloted LESS trajectories in which 
the automatic LESS guidance -parameter values listed in table I were used. 


20 








LESS EVALUATION 


A pilot control station mounted on load cells for transmission of kinesthetic control 
inputs appeared to provide a valid simulation of a simplified lunar escape vehicle. During 
the checkout operations the test subjects were able to adapt to the equipment with a mini- 
mum of instruction and training. Alterations and additions to the control station hardware, 
to the attitude controller characteristics, and to the display instruments were made during 
the progress of the evaluation. The changes and reasons for making the changes are dis- 
cussed in the following sections. 


Test Environment 

Several tests were made in which the pilot was in a pressure suit with a simulated 
portable life support system (PLSS). This configuration is illustrated in figure 13. The 
earth weight of the simulated PLSS , which is the load that a pilot might feel during a 
lunar escape lift-off, is approximately 311 newtons (70 lbf). For this photograph the 
pilot was posed to the left of his proper position in order to permit a view of the display 
panel from the available camera position. 



L-69-5487 

Figure 13.- A view of the pilot during the tests which 
included- a pressurized spacesuit. 
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Two hardware alterations were necessary in order to make the pilot control station 
compatible with pressure-suit operation. A handlebar, as shown in figure 8, was added 
to the handrail. The handlebar permitted the pilot to stand at the appropriate station 
c.g. position and maintain two-hand contact with the handrail. Also the thrust control 
switch could not be positively located by the pilot. The switch box was reoriented and 
the length of the switch lever was increased. The modified switch is the one shown in 
figure 8. 


Controllers 

The electrical dead band used in the axes of the three-axis controller appeared to 
affect control performance by the pilot when the controller was used as a simple on-off 
device. There was a degree of pilot uncertainty as to whether he had exceeded the dead 
band when he attempted to make a short -duration control input. The magnitude of the 
required control motion also caused an increase in the time required for the pilot to ini- 
tiate a desired control input and sense the response. The controller dead-band effect 
was most detrimental when the pilot attempted a small yaw control input. Substitution 
of the single-axis yaw control stick was a definite improvement for the bang-bang type 
of yaw control input which was found to be preferred by the pilots when using kinesthetic 
control for pitch and roll during lift-offs. This control stick was particularly appropriate 
for the tests in which the pilot was in the pressure suit. The mechanism of the single- 
axis controller eliminated a requirement for control "feel" and permitted a more positive 
actuation of the stick-microswitch than was possible with the three-axis controller. 

mation Display 

Several alterations were made to the information display group. The final display 
group is shown in figure 11. The successful completion of a run was stressed rather 
than the fuel quantity used. Therefore the fuel-percentage indicator was eliminated. 

The fuel quantity was recorded, however, and can serve as design information. 

As the tests continued, it became apparent that time was not always a valid situa- 
tion clue to the pilot. If the vehicle had not accelerated to the proper velocity in the 
designated time, then the required orbit would not be achieved. The accumulated thrust- 
axis velocity was presented to the pilot by the large digital display at the top center of 
the panel. Several sizes of digital display numerals were tried before the test subjects 
were consistently able to read the displayed numbers without undue diversion of their 
attention from the 8-ball. The 2.29- by 1.52-cm (0.90- by 0.60-inch) size shown in fig- 
ure 9 was satisfactory. The dial shown to the right provided a differential readout of 


23 



I III! 1 1 II ■ IB III I I I 


II I 


IB 1 1 II 


the velocity (with event markers). When the velocity approached that required for the 
pitch maneuver, the needle on the dial would begin to move. When the needle reached 
the first peripheral marker, the pilot would change the thrust and initiate the desired 
pitchover. The vernier needle was automatically reset at the end of its first traversal 
and did not move again until just prior to thrust-off. The second marker was used to 
indicate velocity for thrust-off. In most cases a single marker would be sufficient for 
all events. This differential readout enabled the pilots to execute the control events 
within 1 to 2 m/sec (3.3 to 6.6 ft/sec) of the scheduled velocity values. This range of 
error was a function of pilot reaction time, vehicle velocity, and vehicle thrust. 

The test subjects tended to lose attitude orientation as they scanned the various 
items on the instrument panel. In particular a requirement to scan upward appeared to 
induce orientation problems. A mirror system was installed (see figs. 8 and 11) so that 
the subject could scan up or down to observe the digital display. All the test subjects 
considered the mirror system a very desirable addition to the panel. 

The marking on the face of the 8 -ball was the source of some difficulty during the 
tests. All of the escape profiles tested began with a pitch attitude of 0° between the 
vehicle thrust axis (Zg) and the local vertical. At this 8-ball position (see fig. 14(a)), 
the top half of the ball is white and the bottom half is black. The line of demarcation is, 
of course, the 0° pitch mark. With the 8-ball at this position, the changes in pitch, roll, 
and yaw angles were readily interpreted and attitude orientation was easily maintained 
(compare figs. 14(a) and 14(b)). At a later point in the escape profiles, the desired pitch 
angle was in the vicinity of -90°. At this angle the pilot's somewhat unfamiliar view of 
the 8-ball would appear as shown in figure 14(c). If oscillating changes in attitude 
occurred simultaneously in two or three axes, as shown in figure 14(d), the pilot could 
lose attitude orientation. In order to avoid this situation, a flip circuit was incorporated 
in the 8-ball reference circuit. As the pilot rotated in pitch from 0° toward the -90° 
attitude position, the 8 -ball could be flipped some predetermined number of degrees so 
that the white-black portion of the 8-ball face could again be used for a reference. The 
flip circuit was activated as the vehicle pitch attitude passed through -60° during the 
90° pitch maneuver. The change at this stage did not disturb the subject and was com- 
pleted well before the pitch maneuver was to be terminated. The advantage gained by 
this procedure can be estimated by comparing figures 14(b) and 14(d). In a real situa- 
tion, it would probably be more appropriate to mark the -90° 8-ball pitch attitude posi- 
tion in a more acceptable fashion. 
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Figure The attitude reference display as seen ty the pilot. 




CONCLUDING REMARKS 


A lunar-escape-system simulator was designed and fabricated to be used as a 
research tool to study problems related to lunar escape and simplified manual guidance 
and control. A pilot station mounted on load cells for transmission of kinesthetic con- 
trol inputs to a digital computer program was shown to provide a valid simulation of a 
simplified lunar escape vehicle. 

During checkout operations the test subjects were able to adapt to the equipment 
with a minimum of training and instruction. 

Several hardware tradeoffs, primarily of the instrument display of situation infor- 
mation, made during the checkout resulted in an improved system concept applicable to 
the design of a prototype lunar escape vehicle. However, the physical size, arrangement, 
and complexity of the information display group were determined to be important factors 
affecting workload and attitude orientation for the lunar escape task. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va., February 19, 1970. 
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APPENDIX A 


LESS AXIS SYSTEMS AND TRANSFORMATIONS 


Body forces and moments acting on the LES are summed in a body-axis system 
( x B’ Y B’ Z b) with origin at the vehicle instantaneous center of gravity (see fig. 15 in appen- 
dix B). Then for the trajectory and orbit calculations, velocity components are trans- 
formed to an inertial system ^Xj,Yj,Zjj and to a local-vertical system (Xlv>Yly,ZlyJ. 
The local-vertical system is used as a reference for a number of simplified LESS guid- 
ance and control schemes. 


Transformation From Body to Inertial System 

The body-axis components of the linear velocity of the LES (u,v,w) are transformed 
to inertial velocity components (xpYpZj) by the matrix operation 


*1 



u 





*1 


fijj 

V 

Z 1 



w 


(i=l,2,3; 3=1 ,2,3) 


(Al) 


where the elements a^ are direction cosines which were generated in the computer 
program (for the (p,ip,6 order of Euler rotations) as follows from the body-axis angular 
velocity components (p,q,r): 


a ll = a 12 r “ a 13 q 
a 21 = a 22 r " a 23 q 
a 31 = a 13 P ~ a ll r 
a 22 = a 23 p " a 21 r 
a 13 = a ll q ~ a 12 P 
a 23 = a 21 q “ a 22 p 


a 31 = a 12 a 23 “ a 22 a 13 
a 32 = a 21 a 13 " a ll a 23 
a 33 = a ll a 22 ' a 21 a 12 

> (A2) 


Further discussion of this generation technique is given in reference 11. Then by refer- 
ring to a standard matrix array of direction cosines (for the order), the Euler 

angles can be extracted as 
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p = sin -1 ^aj2^ 


6 = cos 


-If a ll 


cos p, 


. -If a 32 

cp = -sin 

\cos p. 




(A3) 


These expressions are used in the LESS program to drive an inertial attitude display 
(prototype LM 8-ball having a <p,P,6 gimbal order). Drive signals for an 8-ball which 
is slaved to the local vertical are discussed subsequently in this appendix. 


Transformation From Inertial to Local -Vertical System 

The inertial velocity components Xj, y^, and Zj may be transformed to a 
rotating-local-vertical basis by 


_ *LV 




r -| 

h 

*LV 

= [ C i3 




_ Z LV_ 


_ z i_ 


(A4) 


where the elements Cjj are direction cosines generated as follows by using the angular 
velocities 77 and /3: 


C 11 = ~ C 31 V 

c 13 = c 21 c 32 

" C 22 C 31 

c 12 = " c 32^ 

C 23 = c 12 c 31 

~ C 11 C 32 

C 21 = c 31^ 

c 33 = c ll c 22 

~ C 12 C 21 

c 22 = c 32^ 



C 31 = c ll^ " c 21^ 



^32 = c 12^ " c 22^ 


J 


(A5) 


The angles 77 and /3 are described as the down-range and out-of-plane angles, respec- 
tively, and are defined by 
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7 ] = tan 



J3 = 



(A6) 

(A7) 


where R x ^, Ry^, and 
cussed in appendix B. 



are inertial components of the position vector 


Rq dis- 


Transformation From Body to Local-Vertical System 

By using the preceding two transformations, the body-axis velocities u, v, and w 
can be transformed directly to the local-vertical system by 


— . — 


— — 


— — 

X LV 


u 


u 

^LV 

= Cc.-I Ta-H 

L 1] J LX 

V 

= M 

V 

_ Z LV_ 


w 


w 




(A8) 


where the direction-cosine elements 
bn = c 11 a 11 + c 12 a 21 + c 


b^ have the typical form 
13 a 31 


(A9) 


It is not a simple matter to extract the Euler angles from the b^ elements, so 
the following method is used. The components of local-vertical rotational rate can be 
expressed in the inertial system as 


W LV ,Xj 


£ 

W LV >yi 

■ H 1 

V 

U> T T T 

LV,Zj 
L Li 


0 


(A10) 


where Q indicates a matrix inverse, 
referred to the inertial system by 


i 

B,xj 


p 

B,yj 

■M 

q 



r 

B > Z 1 U 




Similarly, the body rotational rates can be 

(All) 
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APPENDIX A - Concluded 


After summing the two sets of rotational rates in the inertial system, the combined rota- 
tional rates may be referred back to the rotating-local-vertical system by 


r* x 


r 

p 


y 

a 

LV 

COy 

y LV 

•w"‘ 


q 

- M" 

V 

> 

U)rr 

_ Z LV 



r 


_o 

J 


(A12) 


which reduces to 


Wx LV 


— — 
P 




rv-) 

/3 

w yLv 

= 

q 



-1 

V 

co„ 


r 




0 

_ Z LV_ 








(A 13) 


By using these rotational rates, Euler angle rates (indicating body rotational rates with 
respect to the local vertical for the cp,ft,6 order) can be constructed as 


.? 

e 


cos cp' 

y LV 

0)„ sin cp' 

LV 


COS ft' 

cos ft' 

. ' 

ft 

= 

cd v sin cp' + cos cp' 

y LV Z LV 

. » 

<p 


• * # T 

c -6 sin ly 
X LV 


(A 14) 


These expressions can then be integrated to obtain Euler angles which relate the body 
rotations to the local-vertical system. The integrated expressions are used to generate 
drive signals for an 8-ball which is slaved to the local vertical. 
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APPENDIX B 


DEVELOPMENT OF LESS MOTION EQUATIONS 

This section deals with the formulation of motion equations for a body whose mass, 
moments of inertia, and center of gravity vary with time. These variations can occur 
when the propellant mass is a large percentage of the total mass of the body and the rate 
of propellant usage is high. Except for these mass changes, it is assumed that the body 
is essentially rigid. 


Derivation of Force Equations 

The force equations are obtained from the special case of the momentum theorem 
(see ref. 12), which states that the sum of all external forces on a rocket mass is equal 
to the sum of the time rate of change of the total linear momentum and the rate at which 
the exhaust-gas particles are transferring momentum out of the rocket system; that is, 

Z ? ext = - H* + *e) < B1 ) 

where v is the velocity of the rocket mass m and where v e is the exhaust velocity 
of the gas particles with respect to the mass m. After performing the indicated opera- 
tions, this expression reduces to 

^ B ext = ~ mv e (B2) 

In the lunar environment the significant forces acting on a rocket are thrust and the 
gravity. To separate these two forces and determine an expression for the thrust of the 
rocket, the test setup described in the following paragraph is considered. 

A single rocket is held motionless in a test stand such that the thrust axis of the 
rocket is horizontal. During a firing test, v = 0, and the external forces are those due 
to gravity, the pressure at the exit nozzle of the rocket main thruster, and the restraint 
of the test stand. Equation (B2) then becomes 

+ P e A e + f R + - 0 (B3) 

where the first term represents' the gravity force, p g is the directional pressure acting 
normal to the cross-sectional area of the exit nozzle (A e ), and Fr is the total restraining 
force of the test stand. 
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APPENDIX B — Continued 

Next, only the horizontal components of the vectors which make up equation (B3) 
are considered - namely, 

P e A e + mv e = -F h (B4) 

where is the horizontal component of Fr. The rocket thrust F is equal and 
opposite to Fft, and thus the terms on the left in equation (B4) represent the magnitude 
of the rocket-engine force (or gage thrust). Because p g and v e are always parallel 
to F in any orientation, equation (B2) can now be written as 

F + = mv (B5) 

where Rq acts toward the center of the moon; F is defined further in the next section. 

Force Equations Expressed in Body-Axis System 

The origin of the body-axis system (point O) is assumed at the instantaneous c.g. of 
the LES (see fig. 15). The initial c.g. location is designated point G and used for a point 
of reference. The following vector relations are used to define the LES forces in the 
body-axis system (components shown in brackets); 


<1 

o 

II 

x 

(B6a) 

?go = [ Ax > A y> Az ] 

(B6b) 

« = [p,q,r] 

(B6c) 

oU X ?qq = jjq Az - r Ay, r Ax - p Az, p Ay - q AxJ 

(B6d; 

(— r m j = [Ax, Ay, Az"l 

\dt Body L J 

(B6e 

v GO = (| r Go) B0dy + ( S X r Go) 

(B6f 

v O = V G + V GO 

(B6g) 


In equations (B6) , Vq and Vq are the total linear velocities of point O and point G, 
respectively, with respect to the center of the moon, and oJ is the total rotational rate 
of the vehicle about point O; the subscript Body indicates that the differentiation in 
equation (B6e) is performed with respect to the body-axis system. The time derivative 
of Vq is obtained and substituted for v in equation (B5); the resulting vector equation 
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APPENDIX B — Continued 



Figure 15*- LES position vectors with respect to points 0 and G 
and the center of the moon. 


is then resolved into the following body-axis components to give 
for the body forces measured at point G: 

u' + w'q - v'r - Ax(q^ + r%) + Ay(pq - r) 

v + u'r - w'p - Ay(r^ + p2) + Az(rq - p) 

w + v'p - u'q - Az(p2 + q2) + Ax(rp - q) 


F* 


X B 


F* 

y B 

= m 

CQ 1 
. 1 



simplified expressions 


+ Az(pr + q) 
+ Ax(pq + r) 
+ Ay(rq + p) 


(B7) 


The components F^. , F*g, and f£^ are thus body-axis components of the vec- 

tor F (which is a combination of the thrust and the gravity forces). Equation (B7) is 
simplified in that all terms containing time derivatives of Ax, Ay, and Az have been 
omitted as insignificant for the LESS study. 

To relate F Xb , Fyg, anc * *° measures at P oin t O (origin of body-axis 

coordinates), the following relations are used: 
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v Q = [u,v,wj (B8a) 

u = u' + q Az - r Ay (B8b) 

v = v' + r Ax - p Az (B8c) 

w = w' + p Ay - q Ax (B8d) 


The u, v, and w expressions are solved for u', v', and w' , respectively, and also 
differentiated. The resulting expressions are substituted into equation (B7) as indicated 
to obtain the following familiar equation 


F* 

X B 


u + wq - vr 

f ; b 

= m 

v + ur - wp 

jv 


w + vp - uq 


(B9) 


The left side of equation (B9) can be expressed as 


F* 

X B 


YB 


J B 




O/ 


Fy B + b 23 (-^ 




m 


F z b + b 33/ ,J Tl m 

\ R o) 


(BIO) 


where Fx B » Fy^, anC * ^ Z B are k°dy-axis components of the main thrust F and 

the direction-cosine elements b^g, bgg, and bgg resolve the gravity force into the 
body-axis system. No additional thrust terms to represent attitude jets are included 
here because it is assumed that the jets will always be operated to produce force couples 
and thus any asymmetry due to the c.g. shift is discounted. 


The main thrust F is further defined by 


‘V 


f £x b 

F yB 

= 

f *yb 

jV 


-F 


(Bll) 
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APPENDIX B - Continued 

where F is the gage thrust and £ X g and £y^ are thrust misalinement angles (see 
fig. 16). Note that small-angle approximations have been assumed in equation (Bll) 
because and are not expected to exceed ±5° even when the main thruster is 

gimbaled for control purposes. In the LESS studies £ X g and £yg are generated 
either as a pilot’s controller input (when using a gimbaled thruster) or as a fixed error 
(when using a fixed thruster) . 



Figure l6.- Sketch showing LES thrust geometry. 


Velocity-Along-the-Thrust-Axis Equation 

During initial LESS studies it was assumed that a simple integrating accelerom- 
eter would be fixed at point G on the LES to measure only the accumulated velocity in the 
direction of thrust. The following equation was used in the LESS computer program to 
approximate this velocity reading: 




(B12) 
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Moment Equations 


The moment equations are written for body-axis coordinates whose origin is at the 
instantaneous c.g. of the LES (point O). The angular momentum about point O can be 
expressed (see fig. 15) as 


OO f— 

N 0 =I ( R i " R o) X m i( R i ~ R o) 

i=l 


(B13) 


where the vectors - R q) can also be expressed as ^rj - r^) with respect to the fixed 
body point (G). These vectors have the following components in the body-axis system 
(components shown in brackets): 





(B14a) 



(i=l,2,...,°°) 


(B14b) 


?Go = [Ax, Ay, Az] 


(B14c) 


R i 

R i 




? i “ ? GO = [ x i - Ax ’ y i “ Ay ’ z i " Az ] 

(xi - Ai) + q(z t - Az) - r(y* - Ay) 

(yi - Ay) + r(xi - Ax) - p(z t - Az) 

(zi - Az) + p(y t - Ay) - q(xi - Ax) 


(Bl4d) 


(B14e) 


With pilot movements disregarded, the LES is considered a rigid body, so jq, , and Zj 
are zero. Also, Ax, Ay, and Az are assumed zero (see force equations), but Ax, 

Ay, and Az must be retained in the moment equations for control purposes. With these 
assumptions, body-axis components of the angular momentum about point O can be 
expressed (after some algebraic manipulations) as 



OO OO 

- ^ mi(zi - Az)(xi - Ax)r - ^ m^ - Ay)(xi - Ax)q 
i=l i=l 


(B15a) 
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OO 

N y B = I m i[( z i ■ Az ) 2 + ( x i - Ax )" 

i=l 


OO OO 

- ^ m^z* - Az)(y t - Ay)r - ^ m^x* - Ax)(y. - Ay)p (B15b) 

i=l i=l 

OO 

Nz B = I m i|( Xi “ AX ) 2 + ( y l ' Ay ) 2 ] r 
i=l 

OO oo 

- ^ m i (y i - Ay)(zi - Az)q - ^ - Ax^zj - Az)p (B15c) 


i=l 


i=l 


By substituting appropriate inertia symbols for the various summation terms and 
assuming that the LES mass will be distributed symmetrically with respect to the 
X B Z B -, and Y B Z B -planes, the momentum components are reduced to 


n *b 


^■XxP ■ *xz r 


11 

HH 

_V 


*zz r “ ^XzP 
L J 


The following body-torque equations are obtained from the time derivatives of 
equation (B16): 


(B16) 




Q y B 

11 


L 


IxxP + f xxP - I xz(P ( l + r) + (I zz " Iyy)rq 

!y y q + iyyl + *x Z (p 2 “ r ^) + (*XX ~ Izz) r P 

I zz r + I zz r + I xz (rq p) + (lyy - ^xx)Pq 


(B17) 


Terms containing I xz are assumed negligible and thus omitted from equation (B17). 
The pitching acceleration q can be obtained directly from equation (B17) as 


Q y B “ D 2 

q = —7 (B18) 

L yy 

where 

d 2 = iyy^ + ! X z(p 2 “ r2 ) + ^xx “ I Z z) r P (B19) 
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The components Q x and Q z of equation (B17) can be solved simultaneously to yield 

B "B 


and 


i T xz(Qz b - D3) r zz(Qx B - D l) 
P — . . o + , . o 


^xx^-zz 


- (^xz)^ hsx^zz 


- (hz) 


*xz(Qx B “ D l) I xx( < ^z b “ d 3) 


r = 


where 


and where 


*xx*zz _ (*xz) *XX*ZZ (*Xz)^ 
D 1 = *xxP + ( X zz - I yy) qr “ W*! 

^3 = *zz r + (*yy “ *xx)P q + *xz qr 


(B20) 


(B21) 


(B22) 


(B23) 


For the LES studies , Q v , and Q„ are defined as 

X B y B Z B 


Qx b = F Ay - F 5v B ( z Bj t - Az) + Q JjXb + Q P;Xb 

(B24) 

Q y B ' F{ * B ( z B,t «!j,y B * Vy B 

(B25) 

and 


Qz B = Q J,z b 

(B26) 

where F Ax and F Ay are kinesthetic control torques, Q T 

Qj,y B ’ and Qj,z B 


are attitude -jet torques. and are miscellaneous torques (e.g., torques 

^’ y B 

due to uneven propellant drain), and ^£x b ( z b ^ - Az'j and F£y ^z B ^ - Azj are 

thruster-misalinement -torque terms. The angles and £y^ can be either mis- 

alinement angles for a fixed main thruster or gimbal-travel angles for a gimbaled main 
thruster. The distance z D . is the distance the main thruster is displaced (along the 

ij 

Z B -axis) from the initial c.g. location (point G in fig. 16). 
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APPENDIX C 


LESS INERTIA EQUATIONS 

The LES vehicle moments and products of inertia are calculated with respect to 
the body-axis system (Xb,Yb,Zb) with origin at the instantaneous c.g. Because propel- 
lants comprise over one-half the system mass, there is considerable change of inertia 
during an escape-to-orbit run. Thus, inertia-rate effects need to be included in the LES 
moment equations. 


Mass and Center-of- Gravity Relationships 


For convenience the LES mass is divided into the following five parts: mass of 
basic vehicle (m-^), mass of passenger (m2), mass of control pilot (mg), combined mass 
of four spherical propellant tanks when empty (m^), and mass of propellant in the tanks 
(mg). Masses mj, mg, mg, and 014 are themselves constant, but the inertias asso- 
ciated with them are variable because the LES c.g. shifts significantly during an escape 
run. For clarity, the expression "LES c.g." is used herein to refer to the sum of the 
individual-mass c.g.'s (or the instantaneous total c.g.). 


Body-axis coordinates of the following set of five position vectors p m . indicate 
the c.g. locations of the five individual masses with respect to the instantaneous c.g. of 
the LES: 


- [ ix ' iy ’ ( Z G,1 - 4z )] 


(Cla) 

m 2 * [( X G,2 - ax > Ay ’ ( Z G,2 - iz )] 


(Clb) 

m 3 " [(' X G,2 + 5 x B, 3 ' aX )’ ( 6 y B;3 )’ ( Z G,2 - iz )_ 


(Clc) 

m 4 * [( X G,4 - ax )> (y G ,4 - iy > ( Z G,4 - iz )] 


(Cld) 

m 5 = [( X G,4 - ax > ( y G,4 - ay )- ( Z G,4 ' 6 2 B|5 ' iz ) 

(Cle) 


The constants x^ Yq 4* z q 4, mid so forth, are the coordinates of the various indi- 
vidual masses (identified by numeral in subscripts) with respect to the initial location of 
the LES c.g., and Ax, At, and Az represent components of the LES c.g. shift during 

Axj and ^ - Ayj represent the effective horizontal c.g. 

and 6, r are the components of the actual movement 


flight. The terms f6 x 

\ B,3 

shift of the control pilot; 


. Xb , 3 — y Bj3 

of his c.g. from its initial location. The term 


: B,5 


indicates how far the c.g. of the 
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APPENDIX C - Continued 


propellant mass mg shifts downward as propellant is drained from the tanks. It is 
noted that x^, ^ and ^ appear both in equation (Cld) and in equation (Cle); the 
reason is that the initial c.g. locations of and are coincident. Also, Xq 2 

and Zq 2 appear both in equation (Clb) and in equation (Clc) because of the symmetrical 
location of the two masses m 2 and mg. 


Figure 17 is a partial illustration of the mass (and c.g.) relationships; Xq,Yq,Zq 
is a reference system (with origin at the LES initial c.g.) parallel to the primary 
X b ,Y b ,Zb system. Typical mass locations are shown along with just the position vec- 
tor to mg (mass of control pilot, including his suit and PLSS). For this setup, the 
control pilot (facing in x B -direction) has moved forward (®x B ,3) an d to the right (®y B ,3^5 
this movement has effected a corresponding shift Ax and Ay in the total LES c.g. 
Also, the burnoff of propellants has caused the propellant c.g. to drop the distance 6 Z ^ ; 
from the m^ position (i.e., initial ms position) to the ms position shown; this drop 
has caused the upward Az shift in the total c.g. Thus, a change in the vector p m3 
can be caused both by propellant burnoff and the movements of the pilot himself (during 
kinesthetic control runs). 



Figure 17.- LES mass and c.g. relationships. 


40 



APPENDIX C — Continued 


Inertia Equations 

The general forms of the inertia equations are given by 




(C2a) 


(C2b) 


(C2c) 


(C2d) 


In these equations the primes indicate inertias about individual reference axes X, Y, 
and Z (through each c.g.) which are parallel to Xg,Y B ,Z B . Calculation of these 
"prime” inertias is straightforward except possibly for those associated with mg, 
which are treated in the following section. Because of symmetry assumptions, I X y 
and Iy Z are zero. 

Determination of Propellant Inertias 

The "prime" inertias associated with the variable propellant mass mg (with 
respect to a reference system X,Y,Z through the center of a sphere) are herein con- 
structed by a technique of summing elemental disks stacked in the z-direction (see 
fig. 18). Based on this figure, the following relationships, can be written: 


dV = Trx^dz 

(03) 

dm = a dV 

(C4) 

2 2 2 
x = r g - z 

(C5) 

(^Lsk = \ x2dm 

(C6) 


where dV is the elemental disk volume, a is the propellant mass density, and r g is 
the radius of the sphere. The distance x is the radius of any particular disk; the value 
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C *** Continued 


z 

Figure l8.~ Rropellant-taiik geometry showing 
typical elemental disk. 

of x is determined by z, the distance of the disk from the center of the sphere. The 
distance Zj is the vertical distance from the center of sphere to the top disk (viz., the 
propellant (top) surface). (The value of zj is negative when the propellant tank is more 
than one-half full.) The differential moment of inertia of an elemental disk with respect 
to temporary axis X through the center of the sphere is given by 



dI xxg ~ (^Wk + z2dm 


(C7) 


The sum of all the disk-inertia differentials is then given by 

I i X 5 = |i' x2dm+ i' 22,Sm 

= H J ^rj? - z^jdz + no ^ 1 z 2 ^ - z 2 jdz (C8) 

r s r s 


For convenience, let cr be defined as the initial propellant mass (denoted m 5)0 ^ 
divided by the volume of the sphere. By making this substitution and performing the 
indicated integrations, equation (C8) becomes 



2 3 

m c l-r r z. 

5,o 5 s 16 s 1 


l4 

8 r s 80 r 3 


(C9) 


When the propellant tanks are full, z\ = -r s and the equation reduces to the familiar 
expression 
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APPENDIX C - Concluded 


: xx 5 *"5,0(5 "s) 

Similarly, when z^ = 0 (tanks one-half full), the equation becomes 

T » fl 2\ 2 2 

I xx5 m 5,o^5 r s J 5 m 5 s 


(CIO) 


(Cll) 


where m^ is here equal to one -half of the initial mass m^ Q . 

By symmetry 1^^ = lyy^, but does not involve a transfer term for the ele- 

mental masses. However, I zz is similarly obtained by beginning with the relation 


dI zz 5 ~ (^zzJdisk - |( r s " z2 ) dm 


(Cl 2) 


After all necessary operations have been performed, the result is 



= m 



3 

8 r s z l 




(C13) 
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APPENDIX D 


CALCULATION OF ORBITAL PARAMETERS 

The LESS computer program includes logic and equations for predicting the LES 
orbital characteristics based on the burnout conditions in the LES trajectory. In par- 
ticular, the semimajor axis a is determined from 

a = Kq.BO 

2 - R CLB0( v 0,B0) 2 

d 

where the subscript BO denotes values of the state variables at burnout or orbit injec- 
tion (i.e., when thrust is cut off). Next, the radius of pericynthion r p can be deter- 
mined from the following equation: 


r p ~ a 


( R 0,B0) 2 (v 0 ,h) 


BO 


aju 


The altitude of pericynthion h p is then given by 


h p = r p “ r m 


where r m is the radius of the moon. The altitude of apocynthion h a can now be cal- 
culated from the following expression: 


h a = 2a - r p - r m 

The altitudes of pericynthion and apocynthion are the primary variables considered in 
determining whether the LES orbit is safe and/or satisfactory for rendezvous. The LESS 
computer program also contains equations for calculating orbital eccentricity, semi- 
latus rectum, the central angle between injection and pericynthion, and the out-of-plane 
angle (i.e., the central angle between the LES position and the plane of the CSM). The 
equations for these variables are not given herein because this report is not primarily 
concerned with data. 
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